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ABSTRACT 

Binary pulsar systems emit potentially detectable components of gamma 
ray emission due to Comptonization of the optical radiation of the companion 
star by relativistic electrons of the pulsar wind, both before and after termi- 
nation of the wind. The recent optical observations of binary pulsar system 
PSR B1259— 63/LS 2883 revealed radiation properties of the companion star 
which differ significantly from previous measurements. In this paper we study 
the implications of these observations for the interaction rate of the unshocked 
pulsar wind with the stellar photons and the related consequences for fluxes of 
high energy (HE) and very high energy (VHE) gamma rays. We show that the 
signal should be strong enough to be detected with Fermi close to the periastron 
passage, unless the pulsar wind is strongly anisotropic or the Lorentz factor of 
the wind is smaller than 10 3 or larger that 10 5 . The higher luminosity of the 
optical star also has two important implications: (i) attenuation of gamma rays 
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due to photon-photon pair production, and (ii) Compton drag of the unshocked 
wind. While the first effect has an impact on the lightcurve of VHE gamma 
rays, the second effect may significantly decrease the energy available for particle 
acceleration after termination of the wind. 

Subject headings: binaries: close — gamma rays: stars — pulsars: individual 
(PSR B1259-63) 



Introduction 



Three binary systems containing a massive star and a compact object - LS 5039, 
LS I +61 303 and PSR B1259— 63- have been clearly detected in TeV energy band (see 
http : //tevcat .uchicago . edu/| for the updated information). Whil e the nature of the com - 
pact companion in LS 5039 and LS I +61 303 is not yet established (jCasares et al.ll2005aU b: 



Sarty et al.ll201ll ). the detection of the pulsed ra dio emission from PSR B1259— 63 indicates 



the presence of a 47.7 ms pulsar in the system (j Johnston et al.lll992l ). The pulsar orbits a 
luminous star in a very eccentric orbit with the following orbital param eters: eccentricity 



e = 0.87, period P or b = 1237 d, and semi-major axis a<i = 6.9 AU (see iNegueruela et al. 



201 ll . and references therein). The system displays variable broad band nonthermal ra- 



dio, X-ray and TeV g a mma ray emission close to periastron passage (IJohnston et al 



Uchivama et al. 2009 



Chernyakova et al 



2009 




Grove et al. 


1995 





Aharonian et al. 



2005 
20051 . 



20091 ). which currently lacks successful multiwavelength interpretations. Moreover, the Fermi 



LAT observations of periastron passage in December 2010 have shown that in general the 



well 



see e.g. 



Tarn et al. 


2011; 


Abdo et al. 


2011) 



Recently, optical observations with VLT UT2 Kueyen discovered that the optical star 
LS 2883 corresponds to a late O-star an d has a significantly h igher luminosity of = 
2.3 x 10 38 ergs _1 than previously thought ( INegueruela et al.ll201lh . Because of fast rotation 
the star is significantly oblated with equatorial radius of R eq = 9.7R & and the polar radius 
of -Rpoie = 8.1i? . This leads as well to a strong gradient of the star surface temperature 
with T eq = 27500 K and T pr >i e = 34 000 K. The star rotation axis is inclined by ~ 33° 



in respect to the line-of-sight ( INegueruela et al.l 120111 ). The distance to the system is now 
estimated to be 2.3 ± 0.4 kpc. Moreover, the observations favored an orbital inclination 
value of % ~ 25°, which is remarkably smaller than the previously obtained value of ~ 



35 (jJohnston et al.lll994l ). All these new parameters together should have an important 
impact on the multiwavelength properties of this system. 
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The orbital separation distance, pulsar spindown luminosity and the lack of the accretion 
features suggest a realization of the compactified nebula scenario, i.e. the source contains 
two distinct regions: the re lativistic pulsar wind and the terminated flow ( ITavani fe Arons 
19971 ; iBogovalov et al.l 120081 ). The VHE emission is expected to originate in the post termi- 
nation shock region, and a number of models have been propos ed invoking both hadronic 

(Ka wachi et al. 120041 iNeronov fc Chernyakovall2007l ) and leptonic (IKirk et al.lll999l ; iKhangulyan et al 
20071 ) radiation mechanisms. In the framewo rk of the hadronic sce nario, the two humped 
TeV lightcurve obtained with HESS in 2004 (lAharonian et al.l 120051 ) have been interpreted 
as the enhancement of the production rate due to the pulsar passage through the dense 
stellar disc. However, the recent report of a more complicated TeV gamma ray lightcurve 
by HESS with multiple humps and deeps disfavors, to a large extent, the hadronic scenario 
( lAharonian et al.l 120091 ) . In the case of leptonic origin of the emission, there is a number 
of additional effects, which may sig nificantly affect the pro duction rate of the nonthermal 
emission in the post shock region (IKhangulyan et al.ll2007l ). Importantly, in the leptonic 
scenario, one expects a specific radiation component from the unshocked pulsar wind. Al- 
though, the pulsar wind region is not expected to produce VHE emission, a line-like bulk 
Comptonizat ion component from this regio n in HE band is predicted for isolated pulsars like 
Crab pulsar fjBogqyalov Sz Aharonian 200oh and for the binary pulsar system PSR B1259— 63 



faall fc KirkkoodlBall fc Dodd 



2001 



Khangulyan et al.ll2007l ). It is difficult to overestimate 



the importance of observational proof of this radiation component; it offers a unique oppor- 
tunity of detecting a direct signal from the pulsar wind. The most meaningful constraints 
in this regard can be obtained in the 1 to 100 GeV region. Therefore the observations with 
Fermi and AGILE gamma ray telescopes close to the periastron passage are of great in- 
terest; any result (flux upper limit or detections of a signal) can greatly contribute to our 
understanding of the physics of pulsar winds. 

Since the interaction of the pulsar wind with the stellar photon field does not occur in 
the saturation regime, the production rate is very sensitive to the properties of the photon 
field. In this paper we present new calculations of the rad iation signal from the p ulsar wind 
taking into account the new properties of the optical star ( INegueruela et al.ll201ll ) and using 
results of detailed hydrodynamical modeling of the interac tion between the pulsar and stellar 
winds in PSR B1259-63/LS 2883 faogovalov et al.ll2008l ). 



2. Pulsar Wind in Binary Pulsar System 



In the framework of the generally accepted paradigm (IKennel fc Coronitilll9841 ). pulsars 
launch cold ultrarelativitic winds which are terminated due to external pressure. At the 
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termination shock, the wind electrons can be accelerated to multi-TeV energies. The radi- 
ation of these electrons results in a phenomenon called pulsar wind nebula. Since the wind 
is cold, i.e. particles remain at rest in the wind co-moving system, no synchrotron emission 
is expected from the wind before its termination. On the other hand, the comptonization of 
the ultrarelativistic wind by external radiation fields, through which the wind propagates, 
can lead to detectable gamma ray emission. This effect is relatively weak in isolated pulsars, 
and can achieve a reasonable efficiency only in power ful pulsars, provided that the particle 
dominated wind is formed close to the light cylinder (IBogovalov fc Aharonianl 120001 ) . In bi- 
nary systems, the process operates with an enh anced efficiency thanks to the presence of the 
dense radiation field of the optical companion ( IBall &: Kirkll2000l ; iBall fc Doddll200lf ). The 
interaction rate in this channel depends on different parameters characterizing the system: 
(i) luminosity and temperature of the optical star; (ii) orbital separation and inclination; 
(iii) distance to the system; (iv) size of the region occupied by the pulsar wind; (v) pulsar 
wind bulk Lorentz factor. 

Remarkably, the recent optical observations of PSR B1259— 63 have significantly revised 
the parameters (i), (ii) and (iii) in favor of higher temperature and luminosity, smaller 
inclination angle and further system location. Regarding the wind bulk Lorentz factor, it 
remains a highly uncertain parameter. 

The size of the region, where the pulsar win d can propagate depends on the ratio of 
the wind ram pressures, rj (IBogovalov et al.l 120081 ). Since the pulsar spindown luminosity is 
known, one can estimate the ram pressure of the pulsar wind 



J sd 



pw 



Anr 2 c 



where L s ^ and r are the spindown luminosity of the pulsar and distance to the pulsar, re- 
spectively. It should be noted that this relation ignores the possible effects related to the 
anisotropy of the pulsar wind. Although, the level of the anisotrop y may be quite high at large 
dista nces from the pulsar, e.g. in the case of the Crab-like pulsars ( IBogovalov fc Khangoulyan 
20021 ). in this paper we limit our consideration by the case of an isotropic wind. 



To obtain the ram pressure of the stellar wind, one needs detailed information about the 
properties of the optical star, including the mass-loss rate and wind velocity profile, which 
are currently not firmly established. For the given optical star luminosity, the mass-loss 
rate can be estimated at the level of M = 6 x 10 _8 M o yr~ 1 ( Vink et al. 200oh . Accountin g 
for the wind velocity at interaction point V w < Voo = 1350 ± 200 km s -1 (IMcColluml 119931 ) . 



it is possible to estimate the expected value of the ratio of the momentum flux densities, 
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^-parameter, as follows: 



L sd c „ 1n _ 2 / M \ ( K 



-i 



V = = 5 x lO" 2 £ r ^ . . (2) 

' McV w I 6 X lO- 8 M yr- x y V 1350kms / 

We should note, however, that there are several factors which may introduce significant 
uncertainties in the //-parameter. In partic ular, the wind porosity may lead to an overes- 



timation of the mass-loss rate of the star (jOwocki fc Cohen! 120061 ) and consequently to a 



significant underestimate of the rj parameter value. The opposite situation may occur if the 
pulsar wind interacts with the stellar wind close to the star equatorial plane, where a dense 
Keplerian disk is formed. Since the disk is expected to have a significantly higher density 
than the polar wind, and its typical velocity at distance r (i.e. Keplerian velocity) may be 
as high as 

v disk ~ 200 (——) kms" 1 , (3) 



.10 13 cm, 

the disk effective ram pressure may significantly exceed the polar wind one, i.e. the 7]- 
parameter may be remarkably smaller than the estimate of Eq.(T5]). Moreover, because of 
the disk rotation and pulsar orbital velocity, the structure of the wind termination shock, 
in respect to the observer direction, may be rather different for two pulsar-disk interaction 
points. Because of these uncertainties related to the value of the ^-parameter, below we will 
consider a fairly broad range of the r] parameter. 

In Figure [1] the shapes of the termination shock for three different values of the r] 
parameter are shown. The poi n ts in the figure are from the results of numerical modeling 



performed by iBogovalov et al.l (120081 ) . for rj = 1 (squares), rj = 0.05 (filled circles) and 
T) — 1.1 x 10 -3 (open circles). Here the value of r\ — 1 roughly corresponds to the case of 
the interaction with the clumpy polar wind; rj = 0.05 to the case of collision with the stellar 
wind; and rj = 1.1 x 10~ 3 is a lower limit value, which can be realized if e.g. pulsar wind 
is significantly anisotropic at binary system scales; or if the stellar disk plays an important 
role in the interaction. To simplify the calculations we have approximated the termination 
shock by the following analytical expressions: 
for rj = 1 



r 



3.02 yj{z + 0.425649)(z + 0.92) , (4) 



for 7] = 0.05 
for rj = 1.1 x 10" 



0.33^/(z + 0.167)(^ + 5.25), (5) 



r = 0.75y/(z + 0.0289564)(0.367 - z) . (6) 

Here z is the coordinate along the axis joining the pulsar and the star (it is assumed that 
the pulsar is located at the point "0", and the optical star is located at (z = —1, r = 0)), 



-6 - 



and r is the corresponding cylindrical radius (both coordinates are measured in the pulsar- 
star separation units, d p _ s ). Eq.(j6]) shows non-smooth behavior at z ~ 0.367. It does 
not have a physical meaning, but is the result of the simplification of the procedure which 
negl ects the complex stru cture of the termination shock wave close to the symmetry axis 



sec 



Bogovalov et al.ll2008l for details). 



Since the particles in the pulsar wind move radially, the emission towards the observer 
is produced by electrons which were moving in this direction initially. The expected gamma 
ray flux is calculated as 



dN 



termination shock 



dE^dSdt 



And 2 



dl / de 



-ph 



dEjl 



_ T da dN e 

cos 0) e 



dK 



ph 



dE-v dE P dl de ph dV 



pulsar 



(7) 



Here d is the dist ance to the system; da/dE^ is the differential anisotropi c inverse Comp- 
ton cross-section (lAharonian fc Atoyanlll98ll ; iBogovalov & Aharonianl|2000l ); r is the energy 
dependent optical depth due to gamma-gamma attenuation from the gamma ray creation 
point to the observer; and dN ph / de^dV is the target photon density at the given location. 
The term representing the electron density in the cold pulsar wind has the following form: 



diV e 



dE e dl T mc 



■6\E C - T mc 2 - / dl'E/c 



where T and E are the initial wind bulk Lorentz factor and the electron energy loss rate. 

The integration of Eq.flTJ) is performed over the line of sight from the pulsar location to 
the pulsar wind termination shock. Obviously, the integration path depends strongly on the 
orbital phase. In Figure [U the lines of sight for three different orbital phases (—6, and 6 
day from periastron passage) are shown by dashed lines. 

Another effect, which may lead to an additional orbital phase dependence, is the shape 
of the optical star and temperature change between different regions of the star. To study 
this effect we performed calculation for the precise properties of the star, i.e. assuming 
the star to be an oblate spheroid with a linear gradient of the surface temperature as a 
function of the zenith angel. The orientation of the star in respect to the observer is defined 
by inclination angle, i.e. angle betwee n the star rotation axis and line-of-sight, which was 
assumed to be = 33°, as inferred by iNegueruela et al.l (120 111 ). To fix the star orientation 
an additional angle is required, namely the angle which describes the turn in the plane of 
the sky. This angle was assumed to be a free parameter, and its influence was studied. 
Numerical calculations show (see Figures |2] and [3]) that independently on this parameter, 
the emission is well described by a model with a spherical star of radius R* = 6.2 x 10 11 cm 
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and surface temperature T» = 3 x 10 4 K. Given the uncertainties related to the orientation 
of the star, in what follows we perform calculations for the spherical star with the inferred 
parameters. 

In Figure [2] we show the spectral energy distributions (SEDs) expected at the orbital 
phase corresponding to the periastron passage for 77 = 1 (solid lines); 77 = 0.05 (dotted 
lines); 77 = 1.1 x 10~ 3 (dashed line). In calculations we have adapted the following values 
of the initial wind Lorentz factors: T = 10 4 ,4.6 x 10 4 ,2.2 x 10 5 and 10 6 . In Figure E] we 
show a similar plot, but for the orbital phase corresponding to the epoch of 30 days before 
the periastron passage. For a rather broad range of the wind bulk Lorentz factors around 
r = 10 4 , the obtained flux level is above the Fermi sensitivity level (unless the //-parameter 
is small 77 <C 0.05). 

Due to the orbital phase dependence of several key parameters, in particular, separation 
distance, the location of the termination shock, the gamma-gamma optical depth, and the 
electron - target photon interaction angle, the pulsar wind signal has an orbital phase- 
dependence. Since the lightcurves for different energies are quite similar, we show just a few 
examples. In particular, the lightcurve for 10 GeV gamma rays is shown in Figure HI Here 
we assumed for the initial pulsar wind Lorentz factor To = 4.6 x 10 4 . In calculations we 
used two different values of the //-parameter: 77 = 1 (solid line) and 77 = 0.05 (dotted line). 
We note that this parameter may affect not only the flux level, but also the location of the 
light-curve maximum: in the case of the small value the maximum is located close to the 
periastron passage, while in the case of larger values of 77 the maximum is located a few days 
before periastron passage. 

In the case of large bulk Lorentz factor, Tq ~ 10 6 ~ 7 , the inverse Compton (IC) signal 
from the pulsar wind may appear at energies beyond the range of Fermi/LAT. In this energy 
band, the atmospheric Cherenkov telescope arrays are more appropriate tools for probing 
the wind's Lorentz factor (note that for this specific source, currently only the HESS array 
is able to monitor PSR B1259— 63). In Figure H] we show a light curve of 0.4 TeV gamma 
rays calculated for the //-parameter 77 = 1 and for the bulk Lorentz factor of r = 10 6 . 



3. Impact of the higher luminosity of the optical star 

In addition to gamma-radiation of the unshocked pulsar wind, we expec t gamma rays (at 
highe r energies) from the Compton scattering of shock-accelerated electrons (jKennel &: Coronfti 



19841 ) . If the optical radiation density exceeds the density of the magnetic field, the IC gamma 



ray production proceeds in the saturation regime, thus the increased luminosity of the optical 
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star does not lead to amplification of the VHE g amma ray signal. On th e other hand, for the 
recently reported luminosity of the optical star ( INegueruela et al.ll201ll ). the gamma-gamma 
opacity for VHE photons traveling from the pulsar to the observer, may be as large as 0.5, 
thus emission may be attenuated by a factor of 1.6. In Figure we show the corresponding 
optical depth as a function of the phase for four different energies of gamma rays: E 1 = 0.05, 
0.15, 0.4, 1 and 5 TeV. Note that the gamma ray absorption is strongest at the energy of 
0.4 TeV, while at energies below 100 GeV and at multi-TeV energies it becomes negligible. 
We note however that the actual absorption level depends on the production region location, 
while the calculations in Figure |5] assume that the production occurs in the pulsar location. 
In particular, in Figure H] two lightcurves for 0.4 TeV gamma rays are shown. The thick 
dash-dotted line corresponds to the flux level corrected for the gamma-gamma attenuation, 
while thin dash-dotted line shows the intrinsic flux level. It can be seen that the attenuation 
is somewhat weaker than it is expected from the opacity shown in Figure EJ The reason for 
that is rather simple, namely since the gamma ray production occurs along the line of sight 
in the pulsar wind zone (see Eq.flT])), some gamma ray photons suffer a weaker attenuation 
than the one shown in Figure |5j 

Another important implication of the enhanced optical luminosity of the companion 
star for the production of VHE gamma rays is related to the so-called effect of Compton 
drag which leads t o the reduction of the Lorentz factor of the wind before it terminates. 
Khangulyan et al.l ( 120071 ) have argued that this effect could lead to a deficit of power available 
for acceleration of VHE particle at the termination shock. It has been shown, in particular, 



that if the optical star would have a luminosit y close to 4 x 10 



erg s 



then the tendency of 

reduction of the TeV flux observed by HESS ( Aharonian et al.ll2005l ) close to the periastron 
passage may be explained by the reduction of the ove rall energy of the pulsar w ind due to the 



Compton drag. The stellar l uminosity reported by INegueruela et al. 



(120111 ) is remarkably 



close to the one speculated in Khangulyan et al.l ( 120071 ). thus this effect now becomes more 
relevant to the observed TeV lightcurve. To describe the effect quantitatively, we have 
calculated the energy fraction lost by electrons emitted within the solid angle of it steradian 
towards the optical star. The result are shown in Figure |6] for different values of the initial 
wind Lorentz factor: T = 10 4 (dash-dotted line), 4.6 x 10 4 (solid line), 2.2 x 10 5 (dashed line) 
and 10 6 (dotted lines). The //-parameter was assumed to be rj — 1. The calculations show 
that this effect may lead to a significant decrease of the pulsar wind energy transported to 
the termination shock in the case of high values of the ^-parameter, i.e. r] > 0.1. This should 
lead to a proportional decrease of the VHE gamma ray production. Obviously, the Compton 
drag can be important only if the GeV flux from the pulsar wind is high. This radiation 
should be clearly detected by Fermi/LAT, unless the pulsar wind is strongly anisotropic. 
The detection of the pulsar wind signal with Fermi could allow a quantitative estimate of 
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the expected decrease of the VHE gamma ray production. 

The Compton drag should lead as well to a decrease of the pulsar wind ram pressure in 
the collision region. However, given a relatively weak dependence of the termination shock 
shape on the ^-parameter, it is unlikely that the Compton drag would change significantly 
the pulsar - optical star interaction regime. 



Conclusion 



Motivated b y the recent revision of optical properties of LS 2883, the companion star of 
PSR B1259— 63 ( iNegueruela et al.ll201ll ). we present new calculations of hig h energy gamma 



ray fl uxes using a realistic termination shock geometry as described in iBogovalov et al. 
( 120081 ). Calculations show that the higher optical star luminosity is compensated, to a 



large extent, by the new estimate of the distance to the source. Thus, the new calculations 
of gamma ray fluxes are quite close to previous predictions based on the old gamma ray 
parameters of the optical star (see in iKhangulyan et al.l 120071 ) . According to Figure [21 the 
0.1-10 GeV gamma ray fluxes calculated for epoch s close to the perias tron passage are below 
the current upper limits obtained with EGRET ( Tavani et al.lll996l ). and above the mini- 
mum fluxes detectable by Fermi/LAT for observation time of about 1 month. The pulsar 
wind radiation component can be identified by its distinct spectral shape. Another impor- 
tant feature of this emission is the expected modulation with the pulsar period on the top 
of the smooth orbital phase dependence. Indeed, since the emitting electrons move towards 
the observer almost with the speed of light, the gamma ray signal should have the second 
modulation reflecting the time structure of the striped pulsar wind. Although, a detailed 
shape of the fine lightcurve can be hardly obtained given a lack of any consistent description 
of the striped pulsar wind, a detailed search for this effect loo ks quite importan t for a consis- 



tent in terpretation of the results obtained with Fermi/LAT iTam et al.l ( 1201 ll ); lAbdo et al. 

(boiih . 



The most important implication of detection of this component of gamma radiation 
would be the unique opportunity to measure the Lorentz factor of the pulsar wind. On the 
other hand, in the case of failure of detection of this component at GeV and/or TeV energies, 
the conclusions could be equally interesting and important. The possible reasons of non- 
detection of gamma rays from the unshocked wind could be: (i) an extremely powerful stellar 
wind, i.e. very low values of the rj parameter; (ii) unconventional, i.e. very low (Tq <C 10 ) 
or very large (To > 10 6 ) values of the pulsar wind bulk Lorentz factor. The first condition 
requires the pulsar to interact with the stellar disc all over the orbit. This implies a very 
specific realization in the sense of orientation of the stellar disc (namely the orbital plane 
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and the disc plane should almost coincide), which contradic t s to the current e xp ect at ions 
( IMelatos et al.l Il995l ; iBogomazovl 120051 ; iBogovalov et al.l 120081 ; iKerschhaggll 120101 ) . However, 
we should note that one cannot exclude that the pulsar wind is strongly anisotropic. If so, the 
gamma ray signal should be anisotropic as well. This can be another reason for reduction of 
the gamma ray flux, which unfortunately would make the conclusions concerning the range 
of parameters To and 77 less robust. Finally, one should mention that if the pulsar wind is not 
absolutely cold, electrons in the frame of the wind might have a rather broader distribution. 
This would make the gamma ray spectrum less distinct and smoother compared to the ones 
shown in Figures [2] and [3j 

Regarding VHE energy gamma rays prod uced after termination of the wind, the new 
optical observations of lNegueruela et al.l (120 111 ) imply a significant reduction of the flux of IC 
gamma rays produced by shock- accelerated electrons. All three main factors related to (1) 
the larger distance to the source, (2) the gamma-gamma attenuation, and (3) the Compton 
drag of the pulsar wind work in the same (negative) direction reducing the gamma ray flux by 
a factor of up to 10. Given that the previous studies based on the old optical observations 
already have required a significant fraction of the spin-down luminosity (5%-10%) to be 
released in TeV gamma rays, the revised energy requirements become almost unbearable. A 
possible solu tion to the energy budge t crisis could be the Doppler boosting of radiation as 
suggested in iKhangulyan et al.l (120081 ) . This important issues will be discussed elsewhere. 
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Fig. 1. — The geometry of interaction of pulsar and optical star winds: the pulsar is 
assumed to be located at the point with coordinates r = 0, z = 0, the optical star at 
r = 0, z/d p _ s = — 1. The shapes of t he termination shocks as obtained through the numeri- 
cal modeling (jBogovalov et al.ll2008l ) are shown for the following values of the 77-parameter: 
1] = 1 (squares), 77 = 0.05 (filled circles) and 77 = 1.1 x 10~ 3 (open circles). The analytical 
approximations Eq.()3K]) are shown with solid lines. The directions towards the observer are 
shown with dotted lines for three different epochs: —6, and 6 days to periastron passage. 
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Fig. 2. — Spectral energy distributions of IC radiation from the unshocked pulsar wind for 
the epoch of periastron passage. The calculations were performed for different values of the 
//-parameter: t] = 1 (solid lines), rj = 0.05 (dotted lines) and rj — 1.1 x 10~ 3 (dashed lines). 
The initial pulsar wind bulk Lorentz factor was assumed to be Tq = 10 4 , 4.6 x 10 4 , 2.2 x 10 5 
and 10 6 . The thick dashed line roughly corresponds to the expected Fermi sensitivity for 
O.lyr observation. The lines show model calculations performed for a spherical star with 
radius i?* = 6.2 x 10 11 cm and surface temperature T» = 3 x 10 4 K; the regions filled with 
gray correspond to calculations with oblate star for possible star orientations. 
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Fig. 3. — The same as in Figj2]but for the epoch of 30 days before periastron passage. 
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Fig. 4. — The solid and dotted lines show light-curves of lOGeV emission for two different 
values of the ^-parameter: rj = 1 (solid lines), rj = 0.05 (dotted lines) for the initial pulsar 
wind bulk Lorentz factor of Tq = 4.6 x 10 4 . Light curve of 0.4TeV gamma rays, calculated 
for rj = 1 and the initial bulk Lorentz factor of T = 10 6 , is shown with thick dash-dotted line 
(accounting for gamma-gamma absorption) and with thin dash-dotted line without gamma- 
gamma attenuation. The calculations are performed for a spherical star with radius -R* = 
6.2 x 10 11 cm and surface temperature T« = 3 x 10 4 K. 
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Fig. 5. — The optical depth for gamma-gamma attenuation from the location of the pulsar 
to the observer for different energies of the gamma rays. The lines show model calculations 
performed for a spherical star with radius R* = 6.2 x 10 11 cm and surface temperature 
T» = 3 x 10 4 K; the filled regions correspond to calculations with oblate star for possible star 
orientations. 
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Fig. 6. — The averaged ratio of the electron Lorentz factor F at the termination shock to the 
initial value T . The averaging is performed for electrons propagating within the angle of 60° 
towards the optical star. The ratio is calculated for the different values of the initial pulsar 
wind bulk Lorentz factor: r = 10 4 (dash-dotted line), 4.6 x 10 4 (solid line), 2.2 x 10 5 (dashed 
line) and 10 6 (dotted line). The ^-parameter was assumed to be rj = 1, and the calculations 
are performed for a spherical star with radius R* = 6.2 x 10 11 cm and surface temperature 
% = 3 x 10 4 K. 



